Mean field effects in a trapped classical gas 
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In this article, we investigate mean field effects for a 
bosonic gas harmonically trapped above the transition tem- 
perature in the collisionless regime. We point out that those 
effects can play also a role in low dimensional system. Our 
treatment relies on the Boltzmann equation with the inclu- 
sion of the mean field term. The equilibrium state is first 
discussed. The dispersion relation for collective oscillations 
(monopole, quadrupole, dipole modes) is then derived. In 
particular, our treatment gives the frequency of the monopole 
mode in an isotropic and harmonic trap in the presence of 
mean field in all dimensions. 



The dynamics of Bose-Einstein condensates (BEC) of 
dilute atomic gases are described by the Gross-Pitaevskii 
equation ( jjj and references therein). The main feature 
of this equation is the mean-field term arising from in- 
teraction between particles. Most of BEC experiments 
are in the so-called Thomas-Fermi regime for which the 
interaction energy dominates the kinetic energy, result- 
ing in an inverted parabola density shape of the conden- 
sate. However, the mean-field term is not found only in a 
bosonic gas well below its critical temperature T c . Such 
a contribution from collisions also exists above T c , and is 
even magnified by a kind of Hanbury-Brown and Twiss 
factor. Up to now most BEC experiments have been per- 
formed in the collisionless regime (when the mean free 
collision rate is small relative to the trap frequency) and 
with a negligible contribution of the mean-field of the 
non-condensed atoms. 

In order to specify the role of dimensionality, we 
introduce the adimensional parameter £ = gn/ksT 
i.e. the ratio between the mean field energy and the 
thermal energy. One readily establishes that (30 ~ 
(na 3 ) 1 /3( n A3 B ) 2 / 3 , where \ dB = h(2Trmk B T)- 1 / 2 is the 
de Broglie wavelength and a the s-wave scattering length. 
Consequently for a dilute bosonic gas above the critical 
temperature: <C 1 and results presented in this paper 
are valid as corrections. In the weakly interacting limit, 
the 2D quantity &d is only logarithmically small with 
respect to \ 2 dB and mean field energy can be comparable 
to thermal energy above the quantum transition temper- 
ature (Kostherlitz-Thouless The quantity C,\d is of 
the order of (n\dB) 2 (nl c )~ 2 where the correlation length 
l c = h/^/mgn. Classical description and mean field can 
be used up to the regime where XdB ~ 1/n ~ l c so up 
to Ci-D ~ 1- In the regime where l c is much smaller than 
the mean interparticle separation the gas acquires Fermi 



properties and is called a gas of impenetrable bosons or 
Tonkgas|]. 

Finally, by reduction of the dimensionality the role of 
the mean field even above the critical temperature of a 
quantum transition can be important. Experiments per- 
formed on microchip offer the possibility to investigate 
low dimensional regime as well as experience with 
dipolar trap or/and magnetic trap Q. 

In this paper, our aim is to extend the traditional 
treatment of the bosonic gas above the critical temper- 
ature by taking into account the effect of particle inter- 
actions. The method consists in including the classical 
mean field term, also known as the Vlasov contribution, 
within the Boltzmann equation. So far, collective oscilla- 
tions of a bosonic gas above the critical temperature have 
been investigated without mean field contribution in the 
hydrodynamic regime in Ref. [0J|], and an interpolation 
formula from the collisionless up to the hydrodynamic 
regime has been proposed in Ref. 

In Sec. I, we briefly recall the general framework that 
describes how the mean field is taken into account in 
the Boltzmann equation. The stationary solution is dis- 
cussed in Sec. II. In Sec. Ill, we derive equations of the 
low energy collective excitation of a Bose gas for posi- 
tive and negative scattering length by means of a scaling 
ansatz. We report an interpolation formula from the col- 
lisionless gas to the interaction dominated thermal gas 
(Vlasov gas) in the absence of dissipation. 



I. FORMULATION 

In traditional BEC experiments, the bosonic gas above 
the critical temperature is well-described by the classical 
Boltzmann equation or the Uhlenbeck-Boltzmann 
equation Jl^ ] if experiments are sufficiently accurate to 
measure deviation from the classical distribution. In the 
following, we present the extension of this equation when 
the mean- field is taken into account. 

We consider an ensemble of harmonically trapped ther- 
mal atoms that evolves according to the Boltzmann- 
Vlasov kinetic equation |ll],[l3| : 

df df dU 3f gdn df 
dt dr dr av m dr dv 

where /(r, v, t) is the single particle phase space distri- 
bution function, U — 53,w?rf /2 the confining poten- 
tial, n = j fd D v the density, where g is the mean field 
strength in D dimensions. In 3D, g — 4irh 2 a/m is the 
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strength of the pseudopotential replacing the true two- 
body potential at low energies , with the s-wave scatter- 
ing length a. I co \\ is the collisional integral that describes 
relaxation processes. Note that I co a = in ID because of 
conserved quantities. The Vlasov term (last term of the 
l.h.s. of (0)) is a Hartree-Fock mean-field term |Q and 
is even magnified for point-like interactions by a factor 
of two with respect to the condensate for the same den- 
sity. Indeed, for non-condensed cloud both the Hartree 
and the Fock terms contribute, whereas only the Hartree 
term contributes for the condensate. 

The kinetic equation (Q) is valid for ksT 3> Hu) where 
u> is the typical trap frequency and for a -C ao where 
a = (h/ (mu)) 1 ' 2 is the oscillator quantum length. One 
has to check that the s-wave approximation is valid. In 
3D it requires that a <C X c ib ■ 



reached only for low dimensional system. For intermedi- 
ate g, Eq. (U) gives the proper interpolation between the 
gaussian and the Thomas-Fermi shape. 

For g < 0, the density distribution is sharpened with 
respect to the free gaussian one. Actually, by increasing 
the number of atoms, the spatial extent of the distribu- 
tion is reduced. If attractive forces overwhelm the kinetic 
energy, the cloud should collapse. One can work out the 
criterium for such a collapse in 3D by means of a gaussian 
ansatz |1| and finds a c = 33aoN~ 1 (ao/XdB) 5 - However 
this result is out of the range of validity of the classical 
approximation. 

III. COLLECTIVE OSCILLATIONS OF A 
COLLISIONLESS GAS 



II. EQUILIBRIUM STATE 

At equilibrium, the Eq. (jl|) reads: 

E° ( d h „2 r d/o 2gdn Q df 
» 1 1 1 Tn fir : fin.- 



0. 



(2) 



By multiplying Eq. 
and velocity, we deduce the average size 
axis 101: 



dvi m dri dvi 
by VjVj and integrating over space 
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0. 



along the j 



(3) 



As expected, repulsive interactions {g > 0) favor a reduc- 
tion of the density from the free particle situation. The 
opposite behavior is obtained in case of attractive inter- 
actions (g < 0). One can extract the shape of the density 
by searching for a factorized solution of (g) of the form: 
/o(r, v) = n v (v)n (r). One finds a Gaussian spherical 
distribution for the velocity. The density distribution is 
a solution of the following equation: 



Klnno + 2gno/m = fj, 
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(4) 



where n = J Vjn v d D v — hsT/m. In two limiting cases, 
the solution has a simple form. For g = 0, we find the 
gaussian shape as expected for an harmonic confinement 
without the Vlasov term. On the contrary, in the limit 
in which the interparticle interactions dominate and are 
repulsive, the shape of the cloud is determined by a bal- 
ance between the harmonic oscillator and interactions 
energy resulting in approximately an inverted parabola 
shape for repulsive interactions. This is the same shape 
as found for a harmonically trapped BEC in the Thomas- 
Fermi regime since in this case the mean-field term 
also dominates. Note that this last result can also be 
shown in the classical hydrodynamic regime || under the 
same conditions. Strictly speaking this situation can be 



In this section, we investigate the collective oscillations 
of a Vlasov gas i. e. in the absence of the dissipative term 
(^coii) but with the mean field contribution. 



A. Scaling ansatz method 

We study those modes by means of the scaling fac- 
tor method ||,[ll|-[l9| in D dimensions. We recall that 
in this method the proper shape of the cloud does not 
enter directly in the equations. This is the reason why 
the solutions are equally valid for a Bose gas just above 
the critical temperature as for a classical gas. We make 
the following ansatz for the non equilibrium distribution 
function: /(r,v,<) = /q(R(£), V(t)) with Ri = rij\ and 
Vi = \vi — Aj7-j. The dependence in t is contained in the 
free parameters A^. By substituing this ansatz into Eq. 
(Eh, we find: 



E 



Vi df 
A? dFL 



2g dn df 



1 1, A, dRi 



dV l 



(5) 



This equation can be combined with Eq. (|2|) taken in the 
phase space point (r = R, v = V) in order to replace the 
last term of (||) by a linear superposition of dfo /dRi and 
dfo j dVi . We finally obtain: 



E 



A? 



\Ri I Aj + Us Xi — 



XiUjXj 



dfo 
dVi 



dfa_ 

dR t 

= 0. 



(6) 



This equation provides the constraints on our ansatz. 
The first average moment of R%Vi, namely 
/ RiVi[Eq.@)d D Rd D V/N, leads to a set of Newton-like 
second order ordinary differential equations: 
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A, 



- A, - -3 + 



1 

A? 



AjlljAj 



= 



(7) 



with £ = g(no)/(g(no) + fc^T). To find the excitation 
frequencies of the modes, one must linearize around the 
equilibrium value: Aj = 1. 



B. Monopole mode 

Consider the case of an isotropic harmonic confinement 
(uji = luq). In this case the scaling ansatz is in the kernel 
of the collision integral and provides a solution of (|l|) 
valid also in the collisional regime. An exact solution 
of this mode was first reported in |^0| for the classical 
Boltzmann equation without mean field. 

The small amplitude expansion of (0) gives the fre- 
quency of the monopole mode (also called the breath- 
ing mode): Wov4 + £,{D — 2), for all dimensions and in 
presence of the total effects of collisions (mean field term 
and dissipation via I co \\). In 3D this frequency ranges 
from 2lu in the absence of the mean field term up to 
y/Eu)Q when the mean field dominates, in the latter case 
one obtains the same result as expected for Bose-Einstein 
condensate in the Thomas- Fermi limit 25 1 . In 2D we find 
2wo for the monopole mode a result independent of the 
mean field, a special feature of 2D already investigated in 
Ref. [p2| . In this case, the scaling ansatz provides an ex- 
act solution of (|l|). In ID, the monopole frequency ranges 
from 2wo down to y3wo when the mean field dominates. 
Note that this latter result gives exactly the same fre- 
quency as the one of the monopole mode in trapped ions 
p3[ . For ions the force originates from the coulomb in- 
teraction, this long range force is well described by mean 
field and this is probably the reason why we recover the 
same result. 



C. Quadrupolar mode 

In 2D and 3D, Eq. (0) provides the mean field contri- 
bution to quadrupolar collective oscillations for a colli- 
sionless gas. 

In 2D, Eq. (Q) gives two coupled equations for Ai and 
A2, which, after linearization, yield the dispersion rela- 
tion 
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±((4-0 2 (^+^) 



(8) 



For £ = 0, we recover the single particle excitation fre- 
quency of the cloud : to = 2uji for each spatial direction. 
In the limit £ = 1, this relation can be derived from a 
purely hydrodynamic approach by taking into ac- 
count the mean field contribution in the same limit. For- 
mula (|h provides also finite temperature corrections to 



this regime and the proper interpolation in between those 
two limiting cases. 

For a cylindrical 3D harmonic trap (we denote [3 = 
lo z /lo^_) we find the eigenfrequencies of mode M = 
(coupling between quadrupole and monopole modes): 



1 



4 + 4/3 2 - £ 



±([4 + 4/3 2 - £] 2 + 8/3 2 [-8 + 2£ + £ 2 ]) 1/2 



(9) 



and the frequency of the quadrupole mode with az- 
imuthal quantum number M=2: uj 2 /uj\ — 2(2 — £). The 
limit £ ~ 1 gives the formulas derived by Stringari |2q] 
for the low energy excitation spectrum of a BEC in the 
limit for which the energy of interaction predominates 
over the kinetic energy. Low energetical collective exci- 
tations (monopole, quadrupole) spectrum of a BEC are 
rather a proof of mean field dominated physics than a 
direct proof of superfluidity in the Landau sense [Q . As 
already pointed out, the validity of our calculation in 
3D is only perturbative (£ <C 1). For an isotropic 3D 
trap (uji — luq), the oscillation frequency split into the 
monopole mode with a frequency lom — 2o->o(l+£/8) and 
the quadrupole mode with a frequency loq ~ 2loq(1 — £) 
as soon as we take into account the mean field. 



D. Dipolar mode 

The dipolar mode which corresponds to the rigid mo- 
tion of the density profile is not affected by the mean 
field. This can be shown by searching for a solution of 
the form /(r,v,t) = f (R(t), V(t)) with R t = n - th 
and Vi = Vi — f]i. Each component r\i is time dependent. 
Following the same procedure, we readily establish the 
equation of motion for 77 by taking the average value of 
Vi. rji + ujfrji = 0. We recover the fact that Kohn modes 
do not depend on interactions. This last result is nat- 
urally unchanged if we take into account the collisional 
integral contribution. 



E. Discussion 

Note that 7 co n = is strictly speaking only applicable 
for a collisionless gas (a — ► 0) and to the hydrodynamic 
regime (collision rate ^> trap frequencies). In between, 
J co ii is negligible with respect to the Vlasov term if a <C 5, 
where a is a critical value for the scattering length. The 
collision integral is of the same order of magnitude as 
the Vlasov contribution when g = <rv 2 t, where a = Swa 2 
is the elastic cross section, v is a typical thermal ve- 
locity and I a typical size. We take t ~ v/uj and find 
a = 0.03Ads(ArfB/oo) 2 . Just above the critical temper- 
ature a is of the order of the scattering length for the 
experiment of Ref. 01 performed on a microchip. For the 
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metastable helium experiment J26| a is slightly smaller 
than the scattering length. In these experiments even if 
they are not in the Thomas-Fermi regime, one can no 
longer ignore mean field effects. The interpolation pa- 
rameter in Eq. ([?]) is the ratio between the mean field 
and the thermal energy: C = gn/ksT. In many BEC 
experiments, Cmax < 10~ 4 which clearly justifies that it 
is neglected. However, in Ref. this ratio is of the order 
of (max ~ 10% and could be increased by a stronger lon- 
gitudinal confinement. In ^6|, this ratio is of the order 
of Cmax ~ 20±10% as a consequence of the huge value of 
the scattering length and the high density of the sample. 

IV. CONCLUSION 

Mean field effects for a bosonic gas above quantum 
transition play an increasing role as the dimension is re- 
duced. This paper deals with the contribution of the 
mean field to the low energetical collective modes of such 
a gas. We derive, even for a collisional gas, the frequency 
of the monopole mode for an isotropic and harmonic 
confinement. Note that results derived in this article 
hold also for two-components Fermi-system as soon as 
the mean field play a role. 

The mean field contribution could be seen directly on 
the equilibrium shape of the gas in situ. Time-of-flight 
measurement may give relevence to a direct observation 
of mean field contribution J2lj . 

Further experiments on microchip || should allow the 
reduction of dimensionality. Dipolar trap can also be a 
good tool H . Those techniques should help for the obser- 
vation of classical mean field above a quantum transition. 
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